Intrinsic electric polarization in spin-orbit coupled semiconductor heterostructures 
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We present Maxwell equations with source terms for the electromagnetic field interacting with a 
moving electron in a spin-orbit coupled semiconductor heterostructure. We start with the eight-band 
kp model and derive the electric and magnetic polarization vectors using the Gordon-like decom- 
position method. Next, we present the kp effective Lagrangian for the nonparabolic conduction 
band electrons interacting with electromagnetic field in semiconductor heterostructures with abrupt 
interfaces. This Lagrangian gives rise to the Maxwell equations with source terms and boundary 
conditions at heterointerfaces as well as equations for the electron envelope wave function in the 
external electromagnetic field together with appropriate boundary conditions. As an example, we 
consider spin-orbit effects caused by the structure inversion asymmetry for the conduction electron 
states. We compute the intrinsic contribution to the electric polarization of the steady state electron 
gas in asymmetric quantum well in equilibrium and in the spin Hall regime. We argue that this 
contribution, as well as the intrinsic spin Hall current, are not cancelled by the elastic scattering 
processes. 



I. INTRODUCTION 

Spintronic is a rapidly developing and important field 
of condensed matter physics. The research is mainly con- 
centrated on effects of electron spin transport, spin accu- 
mulation and spin manipulation in non-symmetric semi- 
conductor heterostructures with strong spin-orbit cou- 
pling. The early predictions 1 '- 1 ' and recent experimen- 
tal observations of the spin-Hall effect^ 5 - have inspired 
a huge number of theory papers^i^^iii^^^^^ 
This research is primarily concentrated on the spin Hall 
current which is a flux of curriers with opposite spins in 
opposite directions perpendicular to the driving electric 
field. This current can be generated, for a example, due 
to the asymmetric scattering; 1 ^! 3 - due to the diffusion of 
the nonequilibrium spi n 13 i 17 or due to momentum de- 
pendent spin-orbit splitting in the band structure jSiLIS 
The latter effect is usually called the intrinsic effect (as 
it is computed with an equilibrium distribution func- 
tion) and it is characterized by the universal spin Hall 
conductivity*^^ 

Several fundamental questions concerning the spin 
Hall effect have inspired a wide discussion in the liter- 
ature. It concerns the definition of the electron spin 
current ) 11 i 12 i 13 i 14 i 15 and the issue of spin Hall current can- 
cellation in the steady state regime^^^ A good basis 
for treating these issue s 14 ' 15 ! 16 is the relativistic Dirac 
equation for an electron interacting with the electromag- 
netic field. This approach reveals a close relation between 
the spin current and the electric polarization vector. It 
also gives a new contribution to the spin transfer torque 
coming from the interaction between the intrinsic electric 
polarization and the external electric fields 

In the present paper, we describe an electron interact- 
ing with the electromagnetic field and moving in a semi- 
conductor heterostructures with strong spin-orbit cou- 
pling. We start with the eight-band kp Kane model and 



derive the expressions for the electric and magnetic po- 
larization vectors as well as the effective kp Lagrangian 
for the conduction band electrons in a semiconductor het- 
erostructure with abrupt interfaces. Using the least ac- 
tion principle, we derive Maxwell equations with source 
terms and boundary conditions for the electromagnetic 
field at the interfaces as well as equations for the electron 
envelope wave function in external fields together with 
appropriate boundary conditions. As an example, we 
consider the steady state of the electron gas in an asym- 
metric quantum well. We compute the intrinsic spin- 
orbit contribution to the electric polarization in equilib- 
rium, and in the spin Hall regime. We argue that the 
intrinsic electric polarization in the spin Hall regime cor- 
responds to an additional spin-orbit energy in the exter- 
nal electric field, and therefore it can not be cancelled by 
extrinsic contributions. Furthermore, the intrinsic spin 
Hall current does not vanish in the steady state regime 
while the vanishing of the spin torque^ 3 - is maintained 
by additional contributions coming from interaction be- 
tween the equilibrium electric polarization and the ex- 
ternal electric field.— We also predict an existence of the 
intrinsic induced magnetic charge Hall current due to the 
inhomogeneous charge distribution in asymmetric quan- 
tum well. 

The paper is organized as follows: in Section [II] we re- 
view some standard features of the Dirac equation. In 
Sec. IIIII we review the properties of the eight-band kp 
Kane model which includes an external electromagnetic 
field and derive the expressions for the electric and mag- 
netic polarization vectors. In Sec. IIVI we derive the ef- 
fective kp Lagrangian for the conduction band electron 
interacting with external electromagnetic field and apply 
the least action principle to semiconductor heterostruc- 
tures. In Sec. |V]we calculate the intrinsic spin-orbit con- 
tribution to the electric polarization of the steady state 
electron gas in an asymmetric quantum well in equilib- 
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rium and spin Hall regimes as well as the intrinsic induced 
magnetic charge Hall current. In Sec. I VII we discuss the 
results and their relevance for the fundamental issues of 
the proper definition of the spin Hall current, and of its 
non cancellation in the steady state regime. 



II. ELECTROMAGNETIC POLARIZATION 
INDUCED BY DIRAC ELECTRON 

We start with the problem of a relativistic electron 
interacting with the electromagnetic field in the vacuum. 
Maxwell equations for the electric field E and magnetic 
induction B can be written (in SGS units) as: 

1 dE Air 

V-E = Airp, V x B — --- — I J , (1) 

c at c 

1 3D 

V.fl = 0, VxE = --, (2) 

where c is the velocity of light. The charge density p 
and the current density J correspond to a single moving 
electron and satisfy the continuity equation: 

dp 

dt 

For a Dirac electron, p = e^*^ and J = ce(^*a 1 i'), 
where e = — |e| is the free electron charge and \& is a four 
component (bispinor) wave function satisfying the Dirac 
equation 

' eV ) * = (ca ■ tt + me 2 /?)* , (4) 



V- J = 0. 



(3) 



ih- 



dt 



a 
cr 
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(5) 



Here m is the free electron mass, <r = {& x , a y , a z } are the 
Pauli matrices, tt = — i?iV — (e/c)A is the momentum 
operator, and V and A are scalar and vector potentials 
of the electromagnetic fields E and B, respectively: 



E = -(dA/dt)- W, B = VxA. 



(6) 



The Gordon decomposition is a way to split p and J 
into convective and internal parts ^ p = p c + pi, J = 
J c + Ji. The convective parts are given by 

where ^ = ty*f3. The internal parts have the form 

an 

Pi = -V • n , Ji = cV x M + — — , (9) 

at 

where electric and magnetic polarizations II and M are 
given by 

„ e?i - , . 
II = ^(-la)^ , 



M = 



2mc 
eh 
2mc 



<r 
cr 



(10) 



The convective and internal densities are separately con- 
served: dp c ^i/dt + V • J c ,i = 0. The internal density 
Pi and current Ji are the densities of the induced elec- 
tric charge and the induced electric current of a moving 
electron. One can then rewrite the first pair of Maxwell 
equations (P) as 

1 3D Air 

V-D = 4tt Pc , VxH = -— + — J c , (11) 

c at c 

where D = E + AttH is the electric displacement vector 
and H = B — AnM is the magnetic field strength in the 
vacuum. In the general case of dielectric and magnetic 
media, vectors D and H are related to vectors E and B 
via 



D = eE + 4ttII , H = B/ix-AirM. 



(12) 



Here e and p are electric permittivity and magnetic con- 
ductivity of the media, respectively, and vectors II and 
M describe polarizations induced by a moving Dirac elec- 
tron. 

Equation dH/dt = J — J c — cV x M is a direct conse- 
quence of the Dirac equation ^ and its complex conju- 
gate. A similar calculation leads to the following equation 
for the magnetic polarization vector: 



dM 



cV x n . 



(13) 



Here 

Jra 



((tt**q:)S* + *£(c*7r<J<)) - cV x II 

2m 



= 7T- [**P2(W*) + (7T***p 2 )*l 

2m 



(14) 
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is the induced magnetic charge current. It satisfies the 
continuity equation 



dp m 
dt 



V • J m = , 



(15) 



where p m = — V • M is the induced magnetic charge 
density. 

Now Maxwell equations for the displacement vectors 
D and H can be written in the symmetric form: 



V D = Airp c , 

V H = -Awpm , 

„ „ 1 dH Air 

V x D = J m , 

c at c 

_ „ 1 DD Air 

c at c 



(16) 
(17) 

(18) 
(19) 



These equations are similar to those presented in Ref. 1X81 . 
Note that the induced magnetic charge density p m and 
the corresponding current density J m also appear in the 
classical electrodynamics of moving medial 



3 



For a Dirac electron, expressions (| 10[) for II and M are 
exact. Approximate expressions for II and M in weakly 
relativistic limit were recently obtained in Ref. [3 How- 
ever, the Dirac equation or its weakly relativistic limit 
can not be directly applied to the case of semiconductor 



heterostructures with strong spin-orbit interaction^ In 
the next section, we start with the eight -band k ■ p Kane 
model, and we derive expressions for electric and mag- 
netic polarizations induced by the spin-orbit coupling of 
conduction band electrons. 



I 

III. KANE ELECTRON IN EXTERNAL ELECTROMAGNETIC FIELD 

The energy band structure of cubic semiconductors near the center of the first Brillouin zone can be described 
within the eight-band k • p modeh 21 ' 22 In homogeneous bulk semiconductor, the full wave function can be expanded 
as 2 ! 

¥(r)= J2 ^(r)\S)u,+ £ £ ^ va (r)\R a )u, , (20) 

jU=±l/2 [_i=±l/2a=x,z,z 

where U\j2 and U_i/2 are the eigenfunctions of the spin operator S = (h/2)&. \S) is the Bloch function of the 
conduction band edge at the T-point of the Brillouin zone which represents an eigenfunction of internal momentum 
1 = 0. \R X ) = \X), \Ry) = \Y), \R Z ) = \Z) arc Bloch functions of the valence band edge at the T-point of 
the Brillouin zone. Combinations of these functions (\R X ) ± i\Ry))/V% and \R Z ) are eigenfunctions of the internal 
momentum 1=1 with projections on the z axis equal to ±1 and 0, respectively (see Ref. l2ll ). Smooth functions 

(f) are components of the conduction band spinor envelope function ^5> c = ( ^-1/2 j > an d ^^vx (r), W^ty (r), 
^vz (r) are x,y, z components of the valence band spinor envelope vector 

= ( *> ) = { ( ^P 2 ) ' ( ^ ) ' ( Pjj 2 ) } • 



A. Basic equations 



In the bulk, the eight-component envelope function ^>(r) = {^ c (r), ^ v (r)} is a solution of the Schrodinger 



equatio n 23 i 24 i 25 



*>| ( • 1 //Kl " ( «>" I" (21) 




T = 2m 

I , „ A_ *A 



(22) 



-a x 



Here the energy of electron states is measured with respect to the bottom of the conduction band E c = 0, E g = 
E c — E v is the band gap energy, A is the spin-orbit splitting of the valence band, k = — iV is the wave vector, and 
P = —ih{S\p z \Z) /m is the Kane matrix element describing the coupling of the conduction and valence bands. The 
parameter a describes the contribution to the electron effective mass m c which is not related to the interaction with the 
valence band, while the k 2 terms for the valence band are neglected. This is the so called the eight-band Kane model 
with dispersion for electrons onlyj 23 i 24 i 25 This model allows to describe electron states with energies in the conduction 
band. It takes into account spin-orbit effects induced by the interaction with the valence band in the presence of the 
structure inversion asymmetry— Bulk inversion asymmetry terms are not included in the consideration. 
We introduce an expression 

and rewrite the second vector equation of (|21l22p as 

A\ _ iA 
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Taking the cross product with <x, we obtain an equation for it x 

2A\ _ 2iA 



<r x = I E g + - (tx*„ + — *„ . (25) 

Combining (f2"4"| and ([23]) , we express the valence band spinor vector ^f v and the vector product er x ^f v as 

*, = a 1 *i+iC 2 [^x*,], (26) 
<r x *„ = -2iC 2 *; + (Ci - C 2 )[<r x . (27) 

Here the coefficients C\ and C2 are given by 

= 3£ g + 2A = A 

~~ 3E g (E g + A) ' ° 2 - 3£ 9 (^ + A)- lZ8j 

They are related to the electron effective mass m c and the electron effective g-factor g c at the bottom of the conduction 
band via 

77? 

— = a + Epd , g c = g e - 2E p C 2 . (29) 
m c 

Here E p = 2mP 2 is the Kane energy parameter, g e = go + g* , go « 2 is the free electron g factor and g* describes the 
remote band contribution. 

To include the interaction with the electromagnetic field we replace hk = — iTiV with ir = — iTiV — e/cA and 
ihd/dt with ihd/dt — eV in Eqs. (|21l22l23p . and we add the respective Zeeman Kane Hamiltonian -ffzeeman*^ 

tfzecman ( ^ ) = ( 2 ^ [f^ ^ 1 ^ (B(t) vf, ^ + ^fi x ^ ) ' ( 30 ) 

Here (is = \e\h/2mc is the Bohr magneton. In order to simplify our consideration, we assume that the only spin-orbit 
contributions come from the interaction between conduction band and valence band states. Hence, the Hamiltonian 
does not include any Rashba terms related to the remote band contributions. 

In the presence of an additional Zeeman Hamiltonian (|3"0|) the decompositions (|26l27p are not exact. To take into 
account first order corrections coming from ([30]) one has to replace with ^/ + ty B in Eq. (|26I2T|) . Here the 
correction term 

* B m fi B (<rB)(C 1 *i + iC 2 [a x *,]) + H i B B x (Ci*i + iC 2 [a x *,]) (31) 
is by a factor of \i B \B\jE g smaller than \Ev, and it can usually be neglected. 



B. Electron steady state in the stationary electromagnetic field 

We would like to compute the electron steady state with energy £ in the stationary electromagnetic field. To this 
end, we replace ihd^ jdt with {£ — eV)^, where V is the scalar potential: E = — "VV. Then, the valence band 
contribution takes the form 

*„(e) = * £ + *b, (32) 
* £ = -tPCi(e)7T* c + PC 2 (e)[a x tt* c ] , (33) 
*s ~ V B (crB)* e +ifi B B x * e . (34) 

Here e = £ — eV, and the coefficients C\(£ — eV) and C 2 (£ — eV) coincide with coefficients C\ and C 2 given by Eq. 
(f2"5)l after replacing E g by E g +£ — eF. The resulting non-parabolic equation for the conduction band spinor function 
reads: 

7V 2^~(ej ^ + 1 ^ (gc( 4m ge) [ °" * " ] + ^f 9e ^ B ^) * c + lM ' B = e * c ' ( 35 ) 
where again e = £ — eV. The energy dependent electron effective mass m c (e) and ^-factor g c (e) are given by: 

TYl 

— - = a + E p C 1 {e), g c (e) = g e - 2E p C 2 {e) . (36) 
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The contribution of *$>b m Eq. ([34)) [the last term on the left hand side of Eq. (|35)) ] gives corrections to the first two 
terms which are proportional to a small factor \pbB\/ E g , and usually it can be neglected. 

For small energies |e| <C E g , the nonparabolic electron effective mass and electron effective ^-factor can be expanded 
near the bottom of the conduction band 

vn Tfi 

g c {E) = g c + a so e , (37) 



m c (e) m c 



where a p is the mass nonparabolicity parameter and a so is the g-factor nonparabolicity parameter closely related to 
the spin-orbit coupling constant ! 20 ' 26 



E n 



2E n 



El (E g + A 
1 



E 2 g (E g 



E V {C\ + 2Cl) , 



2E V C-2.(2C\ — C2) 



(38) 
(39) 

(40) 



Substituting the expansion of (|37|) into (|35|) and neglecting the contribution of ^b of Eq. ([34| , we arrive at 

Note that \l/ c in Eqs. (|35[) and (|40[) is the original conduction band spinor, and hence the normalization condition 
reads / (\^f c \ 2 + (|*,;| 2 )d 3 r = 1. Expressing via \I> C with the help of Eqs. {32 EH [34j) and keeping only the first 
order terms in \e\/E g and \hbB\/ E g , we obtain the approximate normalization condition for ^S c 



In absence of external fields, this condition can be presented as 

m. r 



^ r \ 2 d 3 r = 



m c (£) 



(41) 



(42) 



C. Electric and magnetic polarizations in the Kane model 

The continuity equation ([3| for the charge density p and the electric current density J in the Kane model can be 
obtained directly from the Schrodinger equation (|2"Tj) . A straightforward calculation leads to the following expressions: 

p = e (|* c | 2 + |*.„| 2 ), (43) 

J = — ((7r* c )** c + tf*7rtf c ) + ieP (***„ - *;* c ) + cV x M , (44) 
Zm 

M = -^^** c -»f° Yl *>*«7-<M*I **«.]■ (45) 



Using the decomposition given by Eqs. (|26ti27|) and neglecting the contribution of *&b in Eq. (|3Tj) one can separate 
convective and internal parts of p and J as p = p c + pt and J = J c + Ji, where internal parts pi and Ji are related 
to electric and magnetic polarization vectors P and M via Eq. (9). For the Kane model, we obtain the following 
approximate expressions for p c , J c , II and M: 

phP 

Pc = e|*d 2 + — (**(V¥) - (V**)* c ) + (46) 
' —^Ji ) + — — (*** + ***„) , 



2 V ^ dt J 2 

ehP - - 

n= — (*:* + *** c ) , * = Ci*„ + iC 2 <x x (47) 
J C = ^_iL_ ( *>* c + ( W * c ).* e ) + ^f:(Vy x CT )* C + (48) 



dt dt dt dt 

M = - l -ti B g c {~eVW c ^c - ™ E ^V*^ - x *„] . (49) 

y=x,y,z 
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It is important to note that the convective charge density p c and current J c enter Maxwell equations |2]) for vectors 
D and H which are related to E and B by expressions (fT!?|) . Vectors II and M in Eq. (fT2")) describe contributions 
of a moving Kane electron into the total electric and magnetic polarizations which are not taken into account by the 
material permittivity tensors e and fi. 

Finally, for an electron steady state with energy |e| in a stationary magnetic field one can apply the decomposition 
given by Eqs. ((32l [33l EH) to Eqs. ((Ml SU ES} . Neglecting the contribution of <5 B in Eq. fat]) and keeping only the 
first order terms in \e\/E g and \/j,bB\/ E gi we arrive at the final expressions for the source terms: 

Pc = e|* c | 2 - (** V 2 * C + V 2 *:* c ) + ^.» B V e (*B)V D , (50) 

n = -£|^V|v]/ c | 2 - ^ (*V x tt* c - (tt*)* x <r* c ) , (51) 

Jc = (*>* c + (tt* c )** c ) + 2|?!***(w x a)*, , (52) 

M = -i MB . 9c (e)*: C r* c . (53) 



IV. LEAST-ACTION PRINCIPLE FOR 
SEMICONDUCTOR HETEROSTRUCTURE IN 
EXTERNAL ELECTROMAGNETIC FIELD 

In the previous section, we found the effective charge 
density, the electric current density and electric and mag- 
netic polarization vectors associated to conduction band 
electrons. The alternative approach is based on the La- 
grangian formalism and on the kp analog of the least 
action principle derived in Ref . H3- This approach is par- 
ticularly efficient in applications to abrupt heterostruc- 
tures. It has two main advantages: (i) the variation of the 
action provides equations of motion together with bound- 
ary conditions at heterointerfaces and (ii) incorporating 
of external electromagnetic field is straightforward. 

The time-independent effective mass Lagrangian den- 
sity for Tq electrons with nonparabolicity is given by* 2 - 7 * 

C{£) = £|* c |2_^!^ V vMrO| 2 + £siA(£XM) 
£ S1A (£) = -^-(g e - g(£))VKW x V* c ] . (55) 

In the presence of external stationary electro-magnetic 
field, the Lagrangian density is 

£cI-EM = £EM + + ^Zeeman ■ (56) 

Here the Lagrangian density £em of the stationary elec- 
tromagnetic field takes into account material permittivity 
tensors e and /i characterizing the material properties in 
absence of moving electrons: 

£EM = T— ( E a Ea(;Ep — B a — Bp) . (57) 

8?r V p, af3 J 

The Lagrangian density £(e) can be obtained from Eqs. 
(|54l55p by replacing £ with e = £ — eV and —ih'V with 
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FIG. 1: Sketch of a planar heterointerface between semicon- 
ductor layers A. E^' B + eV and E^' B +eV are the conduction 
band bottom and the valence band top energies, respectively, 
in the regions A and B. These energies are not defined in the 
boundary region II, while the scalar electromagnetic potential 
V is continuous everywhere. 

7r = —ih'V — e/cA. £zccman corresponds to the conduc- 
tion band Zeeman energy 

£zccma„ = -^. 9 e*:(<rB)* c . (58) 

We consider a planar semiconductor heterostructure 
consisting of two bulk-like regions A and B connected 
by a thin boundary region II around the abrupt het- 
erointerface (see Fig. [1]). The envelope function compo- 
nents $ c are defined only in the bulk-like regions of the 
heterostructure ) 27 ' 28 and they obey the boundary condi- 
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tions at z = — a and z — b. The material parameters m c , 
g c , a, g e , E p and E g may abruptly change from the region 
A to the region B, and they are not defined in the bound- 
ary region II. The scalar potential V and the vector po- 
tential A are continuous throughout the heterostructure 
and do not change inside the thin boundary region II. 
The stationary fields E = — W and B = V x A do 
not have any (5-function components, and they are sub- 
ject to the boundary conditions at z = —a and z = b. 
For the sake of simplicity, we assume that the material 
permittivity tensors e and \i are the same in A and B. 
Then, the Lagrangian densities in the bulk-like regions 
A and B are given by Eqs. (|54II58[) with the material 
parameters of the materials A and B, respectively. Note, 
that the energy E of the electron steady state should be 
replaced with E — E^< B while E + E g should be replaced 
with E - E^' B . Here Ef' B and E^' B denote the ex- 
treme energies of the conduction and valence bands in 
the material A and B, respectively. 

The total action in the heterostructure is given^l by 
$ = h f £ei-EAid 3 r + <5>n- Following the approach 
of Ref. T27I one can show that the contribution >Sn of 
the boundary region in such a model depends only on 
the values of $> c at z = —a and z = b. A variation of 
the action SS = with respect to iff* (with electromag- 
netic potentials V and A assumed to be the constant 
functions of the coordinates) in a standard fashion leads 
to the bulk equation for the electron wave function \t c , 
Eq. (|40|) , together with appropriate boundary conditions 
at the heterointerface. The boundary condition parame- 
ters generally depend on the properties of the boundary 
region II. For the "ideal" interface \a + b\ — > (see Ref. 
|27T ). they can be written as continuity conditions at z = 
of the conduction band spinor function \& c = const and 
of the normal projection v T = (tv) = const (here r is 
the unit vector normal to the interface) of the effective 
velocity vector v: 



1 



m c (e) 



7T* r 



(ffc(e) - ge) 
2m 



[a x tt]* c . 



(59) 



The variation of the action SS = with respect to the 
electromagnetic potentials V and A i: (i = 1,2,3) (with 
the wave function ^ c assumed to be the constant func- 
tions of coordinates) leads to stationary Maxwell equa- 
tions for D = eE + 47I-II and H = B/fj, — 4ttM: 



together with appropriate boundary conditions at the in- 
terface: 



(Dt) = const , (H X t) = const . 



(61) 



Keeping only the first order terms in \e\/E g and 
\/j,BB\/E g , we obtain expressions for the source terms 
p c and J c and for polarization vectors II and M exactly 
in the same way as in Eqs. (|50fl53[) of Sec IIIII Thus, 
the two approaches, the approximation of the eight-band 




V D = iTrp c , VxH = 



4tt 



(60) 



FIG. 2: Sketch of a square asymmetric quantum well with 
infinite potential barriers at z = ±L. The asymmetric dis- 
tribution of the electron density |/n(z)| 2 and of the electric 
polarization TL x (z) are shown schematically. The spin Hall 
current J yz and the induced magnetic charge Hall current J m 
flow in the y direction when the electric field E x is applied. 



model and the least action principle for the nonparabolic 
electrons, produce exactly the same results for the bulk 
semiconductor. In addition, the second approach gives 
boundary conditions at the interface for the envelope 
functions and for the electromagnetic field. 



V. INTRINSIC ELECTRIC POLARIZATION OF THE STEADY STATE ELECTRON GAS IN AN 

ASYMMETRIC QUANTUM WELL 



In this section we consider the spin-orbit contribution to the electric polarization which is created by the in-plane 
motion of electrons in an asymmetric square quantum well. In such a structure, the Rashba-type spin-orbit splitting 
of the electron energy levels can appear in absence of external electromagnetic field (E = and A = 0) due to the 
asymmetry of the interfaces at z = ±L i 26 i 27 i 29 When the interfaces are modeled by infinite potential barriers (both 
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in the conduction and in the valence bands), this asymmetry is reflected by the asymmetric boundary conditions)^ 

(62) 



z=±L 



where k = (k Xl k y ) is the wave vector of the in-plane motion, ai^ = ao/t^, ciq = y l ft 2 /2E p m; t + and t~ are real 

numbers. Here we consider only the case of t ± < which corresponds to the electron energy levels E n > at k = 
(n is the number of the electron subband) . The electron wave function can be written as 



where 



*(»•) = J-^/»(*Mp)> ( 63 ) 
m c {£) 



fn(z) — C + exp(ik n z) + C exp(—ik n z) for|z| < L , (64) 

/»(*)= for|z|>L 



describes the electron quantization, and k n = \j2m c (E n )E n /ft . Constants C are determined by the boundary 

conditions given by Eq. (|62[) together with the normalization condition J^ L \f n \ 2 dz = 1. The asymmetry of the 
boundary conditions ([62]) results in the asymmetry of the electron density distribution |/„(z)| 2 inside the well as 
shown schematically in Fig. 

The function <fi(p) = 4>{x,y) describes the in-plane electron motion and satisfies the equation Hn(£)<fi(p) = £<fi{p) 
with effective nonparabolic Rashba Hamiltonian: 

Hr{£) = h2 f n ±£ } + a S iA{£)W x k] z (65) 
with the effective coupling constant approximated as^- 

T) ^ 711 

a slA (£) = (9e - 9c{£)) (\U-L)\ 2 \.f n (+L)\ 2 ) . (66) 

The second term in the Hamiltonian Hr(E) describes the effective spin-orbit interaction caused by the asymmetry 
of the interfaces. In presence of the external electric field E = (0,0, E z ) the effective coupling constant «sia should 
be replaced by asiA + cur, where cyr ; = (eft 2 /4m)a so E z is the Rashba constant. Note that in this case functions 
f n {z) and the energy levels E n should be calculated taking into account the spatial dependence of the scalar potential 
V(z) = —ezE z . As a result, the asymmetry of |/n(z)| 2 is caused by the asymmetry of the boundary conditions as well 
as by the effect of E z . Although hereafter we will assume E z = 0, all results can be readily generalized for E z ^ 0. 
The eigenfunctions ipx.kip) of the Hamiltonian Hr(E) are given by 

e ikp / x 



= vw 1 ■ (67) 



where k — \k\ — + fc 2 , fc + = k x + ik y , S is the cross-section of the quantum well, and A = ±1 correspond to two 
spin-split spectral branches. Their energies can be found from the equation 

m c (E n ) h 2 k 2 

e«.A,fc = E n r + + Aa S iA(en,A,fe)K. (68) 

m c (e nt x,k) 2m c (€ n: x :k ) 

The spin-orbit interaction shifts the energy minimum of the n-th subband from E n to (E n — E n0 ) [see Fig. GJJa)], 
where the energy shift E n Q can be found from the nonlinear equation 

tp a lih{ E n - E n0 )m c (E n - E n0 ) 

Eno = — 2 ■ (69) 

In what follows we neglect the nonparabolicity of the effective mass and of the coupling constant inside the n-th 
subband and assume m c (e rit x,k) ~ m c {E n ) = m n , asiA(e«,A,fc) ~ a S iA(E n ) = a n . 
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Our aim is to calculate the electric polarization vector: 

n = -^Sr^' 2 _ [K(T x fe * c ~ (M,c) * x °"* c] ■ (70) 

In equilibrium, we have II = (0, 0, IT Z ), and for the electron state characterized by the quantum numbers n, A, k we 
obtain: 

TTTt,A,fc/ \ eft 2 a p m c d\f n {z)\ 2 eh 2 a so m c 2 

The polarization is inhomogeneous in the z direction, and it consists of two contributions. The first contribution 
is due to the nonparabolicity of the electron effective mass and to the asymmetry of the electron charge density 
distribution in the well. It does not depend on the in-plane vector, and it is the same for both spin states of the 
electron. By contrast, the sign of the second contribution is opposite for branches A = 1 and A = — 1. It is related 
to the nonparabolicity of the electron effective g factor. Regardless of this difference, both terms appear due the 
interaction of the conduction band electrons with the valence band states. 

To obtain the full polarization created in the well, one has to integrate over the equilibrium Fermi distribution 
corresponding to the Fermi energy tp: 



n * 9 W =Er4/ n»^ fe (z)d 2 fc = X)n?(*) 

n,\ ^ ' n 



(72) 



At zero temperature T = 0, the integration for each occupied n-th subband should be performed over < k < K±, 
where K± are Fermi momenta for both spectral branches for given n [see Fig. |3] (a)]. K± are determined by 

h 2 K 2 , s 

e F = E n + T a n K± . (73) 

2m„ 

If the Fermi level crosses only the lowest spectral branch A = — 1 of the n-th subband, the integration should be 
carried out over K- < k < K + , where 



a n m n 2m n 

K± = ^ ± J -^g- (e F - E n + E n0 ) . (74) 

For the Fermi energies E\ < ep < E2 — E20 only the first electron subband is filled, and both spectral branches 
are crossed by the Fermi level [see Fig. G2V)]. Then, the integration over k and the sum over A = ±1 give us the 
contribution from the n = 1 subband as 

= (zf - E n + 2E nQ ) « + eana ^ / _ + 4 x _ 

47r m oz 2itTi m \ 3 J 

For the Fermi energies E\ — E w < tp < E\ only the lowest spectral branch of the first electron subband is filled and 
crossed by the Fermi level [see Fig. [3Ja)]. Then, integration over k gives us the contribution from the n — 1 subband 
as 



( = ^ {6F _ En + Eno) ?lipl (76) 

47r m V n oz 



ea so m c 2m n ( 5E nQ \ 2 

When the Fermi energy is increased, more subbands give a contribution into polarization, and the final equilibrium 
polarization can be found as Hl q (z) — J2 n n™( z )- 

Let us now consider the effect of the dc electric field E x in the x direction. We deal with the perturbations 
— —eE x x and — —(eh 2 / 1 Am)a so E x a z k y . The second perturbation is related to the dependence of the 
electron effective gt-factor on E x , and it describes an additional spin-orbit coupling. The perturbation related to the 
dependence of the electron effective mass on E x can be neglected as far as we assume E x to be small and the size of 
the sample in x direction to be large. 

The first order correction to the in-plane wave function caused by is 



(77) 
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The correction to the wave function caused by is smaller by a factor of E n / E g , and it can be neglected here. 
However, we shall later consider the first order correction to the spin-orbit energy which corresponds to H^ 2 \ 

We calculate the intrinsic spin-orbit contribution to the electric polarization n x which is linear in the electric field 
E x as 

n.A 



The result is II^z) = ^-IIJ^), where 



{z) = ym c a E 
lomir 

for E n < €f < E n+ i — £?( n +x)o ! 



and 



m , \ e 2 m c a so E x h j2m n ~ ~~ 2 

n x(^ = 77, J (e F - E n + E n0 ) f n (z) 80 

lbmir m n a n V 7i 

for E n — E n o < E < E n . 

Averaging over z, (Tl x ) z — f^ L H x (z)dz/2L, and introducing an electric susceptibility constant k x , we obtain: 

m^) = %^ = £<m. ( 81 ) 

x n 

The contribution of the n-th subband is given by 

2 

<(e F ) = k = - 32w7rj . (82) 

for £"„ < < -En+l — -E-(n+l)0 ! 

and 



?»/ \ e 2 m c a so % 2m n ~ ~~ 

<(e^) = -^ F \ H: 2~( e F -E n + E n o) 83) 

for E„ - £;„o < €f < S« • 

The physical meaning of the finite intrinsic electric polarization H x can be understood if one considers the spin-orbit 
interaction of the moving electron in the external field E x . This interaction is described by the perturbation 

and the respective spin-orbit energy can be calculated as E so = J d 3 rJ2 n k\ (^'*\H^\ 1 ^' K j . It is easy to see that this 

energy is given by E so = — J d 3 rH x (z)E x = —n x E x fl, where il = 2LS is the sample volume. 

The dependence of the intrinsic electric susceptibility k x {^f) on the Fermi energy is shown in Fig. [3fb) . Remarkably 
the electric susceptibility keeps the constant values K x {ep) — ukq when both subbands A = 1 and A = — 1 of the n-th 
quantum size band E n are crossed by the Fermi level and the next I? n +i band is empty. The value ko of Eq. {82]) 
is independent of the subband number n and of the Fermi energy ep, and it is closely related to the universal spin 
Hall conductivity value a yx = J yz /E x = |e|/87r obtained in Refs. I^l7lfl0l. Here J yz denotes the intrinsic spin Hall 
current corresponding to the flux in the positive (negative) y direction of the electrons with spin parallel (antiparallel) 
to z when the dc electric field E x in the x direction is applied (see Fig. 

The above derivation of the intrinsic contribution to the electric polarization H x (z) in the spin Hall regime allows 
us to predict a new effect of the induced magnetic charge current J m in the y direction when the dc electric field E x 
in x direction is applied (see Fig. [J). Indeed, for the steady state electron gas we have dM/dt = 0, where M is the 
magnetic polarization vector, and according to the Eq. (| 13[) the induced magnetic charge current J m can be defined 
as 

J m = -cV x n . (84) 

Here we assume that the sample is infinite in x and y directions, and we do not consider edge effects. As electric 
polarizations Tl x (z) of (I79I80P are inhomogeneous in the z direction, the y component of the intrinsic induced magnetic 
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charge current is given by J m y{z) = —cdU. x (z)/dz cx d\f{z)\ 2 /dz. Averaging over z, we obtain the mean intrinsic 
induced magnetic charge current: 

(J my ) z = J J my (z)dz/2L = a™E x , (85) 

C(^) = %^ =ECM • (86) 

n 

Here <t^(z) is the average magnetic Hall conductivity summarized over all electron states. The contribution of the 
n-th subband <7™ n is given by 

= - e ^^^i-L)?-\UL)?) (87) 
for E n < e F < E n+ i — -E( n+ i)o , 

and 



= -S^^;V 3 P :i( -- R ' + & ° ,(l/ " ( - i)|3 - |/ " (i)|2 > 

for E n — E n o < eF < E n . 

The dependence of the intrinsic magnetic Hall conductivity cr™(e^) on the Fermi energy is very similar to the 
dependence of the intrinsic electric susceptibility K x {ep) on the Fermi energy which is shown at Fig. [H^b). 
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FIG. 3: Sketch of the electron energy structure in an asym- 
metric quantum well (a), and the dependence of the electric 
susceptibility k x = Xl x /E x on the Fermi energy £f in the spin 
Hall regime (b). 



VI. CONCLUSION 

In conclusion, we have considered the eight-band Kane 
model for the conduction band electrons moving in 
the external electromagnetic field and showed how the 
Gordon-like decomposition can be adapted to this set- 
ting. This approach allowed us to derive the source terms 
for the Maxwell equations and the electric and magnetic 



polarization vectors related to a moving electron. We 
have also derived the effective kp Lagrangian for a non- 
parabolic conduction band and in the presence of the 
external electromagnetic field. In this way, we obtained 
boundary conditions for the envelope function and elec- 
tromagnetic fields at the interfaces. These results give a 
solid basis for the analysis of the spin Hall effect and other 
spintronic effects in semiconductor heterostructures. As 
an example, we have obtained the expression for the elec- 
tric polarization induced by the in-plane motion of the 
nonparabolic electrons in the asymmetric quantum well 
with infinite potential barriers and have calculated its 
dependence on the Fermi energy. We have predicted and 
calculated the intrinsic induced magnetic charge Hall cur- 
rent in the spin-Hall regime. We are going to discuss the 
relevance of our results for the fundamental questions 
concerning the spin Hall effect. 

First, we discuss the definition the spin Hall current. 
At present, three different definitions of the spin current 
operator Jij were suggested in the literature. They are 
the following: (i) the conventional definit ion&£ & 1 1 1 2 1 1 3 
Jij = (H/4)((TjVi+Viijj) (where Vi = \/%{dH /dki) is the 
velocity operator for the Hamiltonian H); (ii) the modi- 
fied definition = (h/2)d{fi(jj)/dt proposed in Ref. [Til 
and (iii) the definition = —ih 2 / (Am c )(<jjdi+di<jj) ob- 
tained in Ref. [3 based on the relativistic approach. We 
note that the definition (iii) ensures the relation between 
the anti-symmetric part of the current A a — Sap-vJp-y 
(here e Q( 3 7 is Levi-Civita anti-symmetric tensor) and the 
electric polarization Hq,^ It is also with this definition 



the universal conductance a v 



|e|/87r remains un- 



changed when the perturbation is taken into ac- 
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count. Finally, with this definition the intrinsic electric 
polarization Tl x is proportional to the intrinsic spin Hall 
current J yz , and the latter can be related to the spin- 
orbit energy E so . This fact also plays an important role 
in the discussion of the cancellation (or non cancellation) 
of the total spin Hall current in the steady state regime. 

The steady state regime is only possible when the scat- 
tering is present in the system. In this case, the to- 
tal spin Hall current J yz includes contributions caused 
by the asymmetric scatteringi2i£, and by the generated 
nonequilibrium spin density S y . 13,17 According to the 
"cancellation theorem" of Ref. ]9|, the spin torque T y is 
proportional to J yz , and it vanishes in the steady state 
regime. However, the spin torque in Ref. was calculated 
without taking into account the perturbation which 
gives an additional contribution to T y . This additional 
contribution l/2E x H e z q to the torque T y was considered 
in Ref. 14, and it can be incorporated in the more general 



form of the "cancellation theorem" presented in Ref. HH . 
The direct relation between J yz on the one hand, and 
and E so on the other hand, implies that J yz does not 
vanish once the contribution to E so corresponding to the 
perturbation is taken into account. If we assume 

that this energy is determined by the intrinsic properties 
of the system, and that it is not affected by the elas- 
tic scattering, we conclude that it is exactly the intrinsic 
part of the spin Hall current which remains non vanish- 
ing in the steady state regime. The extrinsic current is 
cancelled in the way discussed in Ref. 1 1 3L and it does not 
contribute to the total electric polarization n^. 

We are grateful to LA. Merkulov for the discussion of 
the Kane model and to V. L. Korenev for the discussion 
of the spin Hall current "cancellation theorem" . A.V. Ro- 
dina acknowledges the support from the Swiss National 
Science Foundation. 
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